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Abstract 

Ezhcla V.V., Kuyanov Yu.V., Larin V.N.. Sivcr A.S. The Inconstancy of the Fundamental Physical 
Constants: Computational Status: IHEP Preprint 2004-36. - Protvino. 2004. - p. 14. figs. 2, tables 3, 
refs.: 11. 

It is argued that the CODATA recommended values of the fundamental physical constants could 
not be used as the reference data in searching the hypothetical space-time variations of the fundamental 
physical constants. 

It is shown that the CODATA data permanently suffers a loss of self-consistency of the released data 
due to unjustified over-rounding of their estimates. 

The simple estimates of the critical numbers of decimal digits that should be saved in the indepen- 
dently rounded correlation coefficients, the average values and uncertainties to save the self-consistency 
is obtained. 

The set of high level quality requirements to the computerized presentation of the numerical data on 
the jointly measured or estimated physical values are formulated. 

It is argued (once again) that the common standard for presentation of the numerical values of 
correlated quantities in publications and sites is urgently needed. 

AHHOTaiJHH 

Eacejia B.B., KyaHOB IO.B., JlapuH B.H., CuBep A.C. HenocTOflHCTBO dpyHflaMeHTajibHbix cbH3HuecKHx 
nocTOHHHbix: BbiuucjiHTejibHbiii CTaTyc: IIpenpHHT H<J>B9 2004-36. - IIpotbhho, 2004. - 14 c, 2 piic, 
3 Ta6n., 6u6jiHorp.: 11. 

IlpHBe,a;eHbi CBH,a;eTejibCTBa Toro, uto peKOMeH/ryeMbie CODATA 3HaueHHH cbyimaMeHTajibHbix cbH3HuecKHx 
nocTOHHHbix Henpuro^Hbi rjik npoBepKH runoTe3bi o bo3mojkhom pa3JiHHHH 3HaueHuii cbyimaMeHTajibHbix 
nocTOHHHbix b pa3Hbix o6jiacTflx bo BpeMeHH h npocTpaHCTBe. 

IIoKa3aHO, uto ny6jiHKyeMbie CODATA Ta6jiuubi 3HaueHuii xax Ha 6yMa»cHbix HOCHTejiax, Tax 
h b sjieKTpoHHOM BH^e, HcnopueHM HexoppeKTHbiM oxpyrjieHHeM uucjieHHbix npe^CTaBJiemiH cpe^Hux 
3HaueHHH, CTaH^apTHbix OTKJiOHeHHii h KOScbcbHUHeHTOB KoppejiauHH. 

IIpeflCTaBjieHbi npocTbie on,eHKH TOUHOCTeii KoppeKTHoro npeflCTaBjieHHH OKpyrjieHHbix cpe^Hnx 3HaueHuii, 
CTaH,a;apTHbix otkjioh6hhh h KoacbcbHirHeHTOB KoppejifliniH. 9th on,eHKH mojkho Hcnojib30BaTb fljia 

KOHTpOJIH KOppeKTHOCTH H COrjiaCOBaHHOCTH 3HaueHHH cbyHflaMeHTajIbHblX (pH3HHeCKHX nOCTOflHHblX. 

CcbopMyjinpoBaHbi npe^jiojKeHHfl no o6hj,hm Tpe6oBaHH\HM k KauecTBy npe^CTaBJieHHa uhcjiobmx 

flaHHMX O COBMeCTHO H3MepeHHbIX HJIH OUeHeHHblX d)H3HHeCKHX BejiHHHHax: hx cpe^HHx 3HaueHHH, 

CTaH,a;apTHbix otkjioh6hhh h KoacbcbHirHeHTOB KoppejiflinrH b ny6jiHKan,Hflx, cnpaBOUHHKax h Ha caiiTax. 
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1. Motivation 



The possible space and time variations of the fundamental physical constants (FPC) contin- 
uously attract much attention of different investigators since the time when Dirac has invented 
the idea. Following a recent review of J. P. Uzan [1], a general strategy for searches of the 
variability can be outlined as follows: 

• The hypothesis of constancy of the FPC can and must be checked experimentally. 

• It only make sense to consider the variations of dimensionless combinations (ratios) of the 
fundamental constants. 

• If the FPC vary, they most probably vary jointly and slowly. This means that to notice 
FPC variations we should: 

— select several well separated space-time regions; 

— measure/estimate as precise as possible physics observables expressed in terms of the 
FPC, that refer to the same space-time region; 

— compare values of constants in the different space-time regions, but extracted from 
the "space-time region dependent" observables with the same current FPC evaluation 
and adjustment methods. 

Let Vx,i denotes the set of FPC related random variables to be estimated and adjusted 
by the method of least squares (for example) on the experimental data at space-time region 
X. This means that after the successful adjustment we will have in the parametric 1^-space 
the vector of average values (Vx,i) and the corresponding covariance matrix Cov(Ux,i, Ux,j), 
characterizing the interior of the "scatter ellipsoid" centered at the end of the vector of averages 



The same ellipsoid can be represented with the help of correlation matrix Cij(X) = Cor(Ux,i, 
Ux,j) and standard deviations Ux,i = y Cov(Ux,i, Ux,i) of Vx,i- 



^2(Vx,i - (V x ,i)) • [Cov(U x ,i, Uxj)}- 1 • (V x ,j - (V x ,j)) < 1. 



(1) 



E 



Vx,i - (Vx,i) 
U x ,i 



• [Cij(X))] 



_! Vxj - (Vx,j) 
Ux,j 



< 1. 



(2) 
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To see the space-time variability we should see the well separation of the scatter ellipsoids in the 
V-space. Let us say that vector V deviates from the scatter ellipsoid obtained for X space-time 
region by Rx(V, (Vx)) standard deviation if 

J2(Vi - (Vx,*)) • [Cov{V x ,i, Vxj)]- 1 • (Vj - (V x ,j)) = R 2 X (V, (V x )). (3) 

ij 

Then it is easy to see that the scatter ellipsoids obtained for the X and Y regions will be well 
separated if in the whole V space we will have 

R X (V (Vx)) + R Y (V (Vy)) > 2, (4) 

that means that the scatter ellipsoids do not intersect. Hence, to be able to notice the variability 
we should have both: accurately estimated average values and corresponding scatter ellipsoid 
for every space-time region where we estimate the FPC. It is the delicate problem as we will 
show further. 

The only and the best known well elaborated procedures to evaluate and adjust fundamental 
physical constants are implemented at the NIST Physics Laboratory [2]. The set of FPC peri- 
odically adjusted at NIST is recommended by CODATA as the reference source of the FPC for 
scientific applications and technology. In any attempt to notice the space-time variability of the 
FPC one cannot avoid the CODATA recommended values, deemed in the physics community 
as the one of the best known set of FPC adjusted in the space-time region where we are. But 
unfortunately it is impossible. Simply because we never had the set of the recommended FPC 
correct enough for the testing their space-time variability. To show this let us select subsample 
of the dimensionless FPC from the CODATA-2002 recommended set [5], say the set: 



Standard FPC name 


Symbol 


Value (2002) 


Uncertainty 


fine-structure constant 


OL 


7.297 352 568e-3 


0.000 000 024e-3 


electron-muon mass ratio 




4.836 331 67e-3 


0.000 000 13e-3 


electron-proton mass ratio 


m e /m p 


5.446 170 2173e-4 


0.000 000 0025e-4 


electron-deuterium mass ratio 


m e /m d 


2.724 437 1095e-4 


0.000 000 0013e-4 


electron-proton magn. moment ratio 




-658.210 6862 


0.000 0066 


muon-proton magn. moment ratio 


V^Iv-p 


-3.183 345 118 


0.000 000 089 


proton g factor 


g p = 2fj, p /nn 


5.585 694 701 


0.000 000 056 



The corresponding CODATA-2002 correlation matrix is as follows: 



Cor(2002) 


OL 


m e /m M 


m e /m p 


m e /m d 






fjp = 2fi p /n n 


a 


1.000 


-0.247 


0.000 


0.000 


-0.003 


0.230 


-0.002 


m e /m M 


-0.247 


1.000 


0.004 


0.004 


0.008 


-0.934 


0.008 


m e /m p 


0.000 


0.004 


1.000 


0.894 


0.000 


-0.004 


-0.046 


m e /m d 


0.000 


0.004 


0.894 


1.000 


0.000 


0.012 


-0.041 




-0.003 


0.008 


0.000 


0.000 


1.000 


-0.008 


0.999 


VhIVp 


0.230 


-0.934 


-0.004 


0.012 


-0.008 


1.000 


0.350 


fj p = 2^ p /fi n 


-0.002 


0.008 


-0.046 


0.041 


0.999 


0.350 


1.000 



This matrix is non-positive definite matrix (it has one negative eigenvalue = —0.000293338). 
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This means that we have no scatter ellipsoid, the corresponding "scatter region" is unbounded 
and the comparison with any other evaluations is senseless. This confusion might be due to mis- 
prints in the resource database as of 2002, but this is not the case. The same situation with 
non-positive definite correlation matrices is present in all releases of the FPC produced by NIST 
and approved/recommended by COD ATA. Further examples of the wrong subsamples of the 
COD ATA recommended FPC see in the Table 1, where we compare data from the last three 
releases (V.3.0, V.3.2, V.4.0). The other examples presented also in our previous papers [6, 7 
on this subject. 



Table 1. Comparison of the selected CODATA:1986, CODATA:1998, and CODATA:2002 recom- 
mended values for the triads of quantities: averages, uncertainties, correlations. 



CODATA:1986 


Symbol [units] 


Value (uncertainty) X scale 


Correlations 


Elementary charge 
Plank constant 
Electron mass 
l/a(0) 


[C] 

h [Js] 
m e [kg] 
a(O)" 1 


1.602 177 33(49) X 10" 19 
6.626 075 5(40) X 10" 34 
9.109 389 7(54) X 10" 31 
137.035 989 5(61) 


e h m e 


0.997 

0.975 0.989 
-0.226 -0.154 -0.005 


CODATA:1998 


Symbol [units] 


Value (uncertainty) X scale 


Correlations 


Elementary charge 
Plank constant 
Electron mass 
l/c*(0) 


[C] 

h [Js] 
m e [kg] 
a(O)" 1 


1.602 176 462(63) X 10 -19 
6.626 068 76(52) X 10" 34 
9.109 38188(72) X 10" 31 
137.035 999 76(50) 


e h m e 


0.999 

0.990 0.996 
-0.049 -0.002 0.092 


CODATA:2002 


Symbol [units] 


Value (uncertainty) X scale 


Correlations 


Elementary charge 

Plank constant 
Electron mass 
l/c*(0) 


[C] 

h [Js] 
m e [kg] 
a(O)- 1 


1.602 176 53(14) X 10" 19 

6.626 0693(11) X 10" 34 
9.109 3826(16) X 10" 31 
137.035 99911(46) 


e h m e 


1.000 

0.998 0.999 
-0.029 -0.010 0.029 



The eigenvalues of these correlation sub-matrices are as follows: 

COD AT A : 1986 {2.99891, 1.00084, 0.000420779, -0.000172106}; 
COD AT A : 1998 {2.99029, 1.01003, -0.000441572, 0.00012358}; 
COD AT A : 2002 {2.99802, 1.00173, 0.000434393, -0.000183906}. 

Definitely something is wrong with the NIST evaluation/adjustment/presentation procedures. 
We suspect that the origin of these permanent confusions is the unjustified independent round- 
ing of the output interrelated quantities: vector of constant estimates, their standard devia- 
tions(uncertainties) and their correlations. 

Superficial independent rounding may lead to catastrophic changes in the connection of 
averages, standard uncertainties and the scatter ellipsoid: the rounded average values may get 
out of the "etalon" scatter ellipsoid obtained after rounding the correlation matrix. The "scatter 
region" may turn to become hyperboloid. From the other hand any numerical calculation is 
performed with rounding or truncating decimal numbers. 
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To preserve the general properties of the FPC data structure, a special quality assurance 
procedures should be developed and applied. In the next section we collect the high level require- 
ments to the set of FPC needed to guarantee the safe and correct usage of this key informational 
resource. 

2. High level requirements to the set of adjusted FPC 

Let us introduce a few special notations and definitions for different sets of FPC to simplify 
formulation and discussions of the requirements. 

V B or "basic FPC" is the set of constants that participated in the fits to the experimental 
data via observational equations. 

V D or "derived FPC" is the set of constants and units conversion factors that are known 
to be function dependent on basic constants. Symbolically V D = F(V B ) and they are evalu- 
ated on the basis of the V B with the proper propagation of the uncertainties with the sufficient 
accuracy to guarantee positive semidefinitness of the derived covariance martix 1 . 

V or "adjusted FPC" is the V B U V D with cross covariances (correlations) added with 
sufficient accuracy to obtain combined covariance matrix as positive semidefinite matrix. 

yR or "recommended FPC" is the V A but rounded by NIST to be compactly presented 
in their publications and as recommended data for science and technology by COD ATA. 

All data sets V 1 defined above have the same pair of structures: 
V 1 = {Average^ 1 ), Covariance^ 1 )} 

or 

V 1 = {Average^ 1 ), Uncertainty^ 1 ), Correlator^ 1 )}. 

Let us call the internal calculational accuracy of numerical presentation of all components of 
the V B obtained from the adjustment procedures as etalon accuracy. 

2.1. Correctness and Self-consistency 

If the adjustment of the constants belonging to V B is successful then we have positive 
definite covariance (correlation) matrix presented with an etalon accuracy, as well as the vector 
of average values. 

We say that the V D , V are correct if their covariance (correlation) matrices are positive 
semi-definite. In other words, we have sufficient internal calculation accuracy to obtain correct 
results. 



x By definition the covariance (correlation) matrix for the jointly measured or estimated quantities is the positive 
semidefinite matrix, moreover if adjustment is performed by the least squares method the covariance (correlation) 
matrix if presented with the etalon accuracy should be positive definite for the successful adjustment. 
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We say that the V is correct and self-consistent if one of two possibilities is true: 

1) V R = V A or 

2) For any subset v(V R ) C Average(V R ) for which corresponding covariance submatrix 
Cov(v(V R )) is positive definite we have 

[v (y R ) - v(V A )]i • [Cov{v{V A ))]^ . [v(V R ) - v(V A )]j < 1 

or 

[v( yR) - v{V A )]i • [Cov{v{V R ))]-f • [v(V R ) - v{V A )]j < 1. 

These conditions guarantee the self-consistency of the V R , e.g. that the rounded and un- 
rounded scatter ellipsoids are well intersected and unrounded and rounded subvectors belong to 
that intersection. 

2.2. Reliability 

We will say that the next release V R Y is reliable if it is correct, selfconsistent, and if any 
subvector v(V R Y ) with positive definite covariance is ended in the point inside the scatter 
ellipsoid for the corresponding subvector of the previous release. For example, for the 1998 and 
2002 releases these conditions will read 

HV R ) - v(V A )]i • [Cov{v{V A ))]-f ■ [v(V R ) - viV&h < 1. 

The reliability indicator is constructed with an assumption that the relative time variation of the 
fundamental constants during two successive sessions of the adjustments are negligible compared 
with the average relative standard deviation of the constants. 

2.3. Availability 

Next important quality indicator we propose is the availability of all data on FPC (average 
values, uncertainties, correlations) in computer readable forms with as maximal as possible 
completeness and accuracy of numerical data. The importance of the availability is hard to 
overestimate in the era of the Web communications and Web and GRID computations 2 . 

It turns out that NIST and COD ATA, in spite of the nicely organized affiliation web cites 
offer the current and archived data on the FPC in the hopeless obsolete manner, as it will be 
shown in the sections to follow. 

2.4. Traceability 

The traceability in the context of usage the recommended FPC is the access to all input 
experimental and theoretical material used in the adjustment as well as detailed descriptions of 
the used procedures needed to reproduce the adjustment independently in case of any suspicions 
on the misprints in the database, ideological or software bugs. 



2 To taste the importance of the availability requirement we will recommend reader to try to check our calcu- 
lations presented in the motivation section, including the correctness of data extraction from NIST publications 
and site. 
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3. Safety rounding off the correlated quantities 



Here we derive a simple sufficient estimates on the accuracy of a safely independent rounding 
off the average values Vi, uncertainties Ui, correlations Cij obtained in jointly measurement or 
estimation procedures with sufficient etalon accuracy. 

Let (Vi, Ui, Cij), i,j = l,...,n be the aggregate of n jointly measured or estimated 
physical quantities, where numerical parts of Vi, Ui are the real numerical vectors, Ui > 0, Cij 
is the real, symmetric, and positive definite matrix with matrix elements bounded as follows: 

Cu = 1 for all i — 1, . . . ,n and \Ci^j\ < 1.0. 

Suppose that for some reason we need to store and exchange numerical data on this aggregate 
rounded to some accuracy A that is lower than the etalon one. 

Let Rij be the "rounder" matrix, such that if it is added to the matrix Cij, the obtained matrix 
Cfj = Cij + Rij will be real, symmetric, positive definite and all |C^j| < 1 are decimal 
numbers with A digits to the right of the decimal point. 

It is easy to see that matrix Rij should have the following properties: 

Ru = forall i = l,...,n and \Ri&\ < 5.0 X 10 _A_1 . 

Let further c\ < • • • < c n , p x < ■ • • < p n , and cf < ■ • • < be the ordered sets of 
eigenvalues of the matrices Cij, Rij, and C^ correspondingly. Then from the Weil's theorem 
for any j = l,...,nwe have the following inequalities [8], [9]: 

Cl + pi < cf < Q + p n . 

>From the Gershgorin's theorem on the distributions of the eigenvalues of the Hermitian 
matrices [8] it follows that 

Pl >_(„-!). 5 • 10-( A + X ) = • 10~ A 

2 

and hence to have the matrix C^ as positive semi definite matrix it is sufficient to demand 

(n — 1) . r , 

< ci - • 10~ A < cf . 

>From the left inequality we have the final estimate for the threshold accuracy index for 
safely uniform independent rounding of the positive definite correlation matrix Cij with minimal 
eigenvalue c a = \^ in 

( n - 1 

'°sio ( jr^c- 



A> = 



mzn 



(5) 



NOTE. According to the Weil's theorem any uniform rounding the off-diagonal matrix elements 
of the positive semi-definite correlation (covariance) matrix is forbidden. 

Indeed, as rounder matrix is traceless Hermitian matrix, it obliged to have the negative 
minimal eigenvalue. Furthermore from the left inequality of the Weil's theorem statement it 
follows that any rounding could lead to the matrix with negative minimal eigenvalue. 
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Now let us clarify to what accuracy we may round off the Vi and Ui in the decimal presenta- 
tions. Let RY be the such "rounding vector" that the obtained rounded vector V/ 1 = Vi — RY 
is still in the etalon scatter ellipsoid. Then from the condition (2) for the components of the 
rounding vector we will have 

R v RY 

ij ■> 

In the eigenbasis of the etalon correlator Cij the expression (6) can be transformed to 

ij mn 1 rn 3 

where L is a rotation matrix. As we try to find the sufficient condition for rounding vector 
components it is enough to demand the validity of (7) for all correlator eigenvalues replaced 
with minimal one. Then the inequality (7) will become 

< (8) 

Inequality (8) means that we can round components independently only inside the maximal 
hypercube imbeded into scatter ellipsoid: 





(9) 

To obtain the accuracy AY for the i-th component that will be sufficient to guarantee that the 
end of the vector V^ belongs to the interior of the etalon scatter ellipsoid it is sufficient to have 

\RY\[uniti] < 5 • 10-( A i +1) [uniti\. 

>From this bound it follows that to have the rounded vector of average values pointing to the 
interior of the etalon scatter ellipsoid one should save 



Ai > A\ = 



1 / n 

- lo Sio 



2 ll) U-Xg in -(Ui/[uniti]y 

digits to the right of the decimal point. 



(10) 



Now let us turn to the rounding of the uncertainties Ui. It is the common practice to present 
the average values and uncertainties with the same accuracy AY = AV . With this rule let us 
rewrite inequality (9) in the form 



lo gio(^i) > lo Sio 




n 



min 
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Taking into account the equality 3 



Lio gl0 (i/i)j + 1 = py - Ay, 

where Pf is the precision of the Ui we will obtain 



PF > 



1 / n 



2 l0Slu V4.A£ 



mm 



One can see that right part of the inequality does not depend on index i, so we can introduce 
P u which is the same for every i: 

pU _ pU 

The equation (11) give the minimal precision that should not be reduced if we adopt the rule 
that accuracy of the uncertainties should be equal to the accuracy of the average values. 

In summary: we have obtained n + 1 reference numbers and AY defining the levels with 
safety independent rounding off the decimal numerical presentation of the interrelated random 
quantities: average values, their uncertainties, and correlations. 

Having these numbers the strategy for the safety independent rounding can be as follows: 

In self- consistent numerical presentation of interrelated random quantities (Vi, Ui, Cij) 
in decimal real numbers the average values Vi and the uncertainties Ui should have at 
least AY digits to the right of the decimal point and the correlation coefficients Ci^j 
should have at least A^ digits to the right of the decimal point. 



4. Do the CODATA 2002 recommended FPC meet the high level quality 

requirements? 

In this section we present some further evidences of violations of the above high level require- 
ments in the recent releases of the CODATA recommended values of the FPC. 

4.1. Correctness & Selfconsistency 

In motivation section we already presented the evidences that the CODATA data on corre- 
lations are incorrect. Here we present an evidence that the average values of the recommended 
FPC are also questionable, because of over-rounding can easily move them out of the etalon scat- 
ter ellipsoid. To check this the whole adjustment process should be repeated with the "etalon 
accuracy". 

It turned out that we managed to collect enough amount of data from the NIST publications 
to reproduce all steps of the evaluation and adjustment of the basic set of constants [10] only 
for the 1998 release. We had obtained the "correct set of the basic constants" using methods 

3 This equality is valid for real numbers only. For the integer number that treated as the numbers with infinite 
precision it is not valid. 
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described by NIST experts [4] 4 and then calculated the threshold accuracies for the elements of 
the correlation matrix, the averages and the uncertainties. The results are as follows: 

^C,min ~ 7.58 • 10 7 , 

A% = 8 (versus A% ODATA = 3), 
P u = 4 (versus pU,CODATA = 2) 

One can see that the CODATA data suffers the loss of self-consistency of the released data due 
to unjustified over-rounding of their results. 

Having the data on the FPC in the "etalon accuracy" we are able to show that the obtained 
estimates for the threshold rounding indices are indeed close to the real situation and should be 
used as regulators for the correctness of the rounding. To show that the rounding procedure can 
move the end of the vector-of-constants out of the etalon scatter ellipsoid we will use the sample 
of constants that was mentioned in [1] as the candidates to trace the large-scale space-time 
variability of their dimensionless combinations: 

Table 2. Selected basic and derived constants from the IHEP adjustment based on the NIST 1998 input 
data. 



Symbol[units] 


Average value Uncertainty 


h [Js] 
m e [kg] 
m p [kg] 
m n [kg] 
[C] 


6.62606875610000 X 10~ 34 5.2200000 X 10~ 41 
9.10938187491360 X 10~ 31 7.2057063 X 10~ 38 
1.67262158291420 X lO" 27 1.3235274 X 10~ 34 
1.67492715608612 X 10~ 27 1.3253602 X 10 -34 
1.60217646198672 X 10~ 19 6.3181739 X 10~ 27 



The corresponding correlation matrix of their uncertainties in the "etalon accuracy" 5 



Cor 


h 


m e 


m p 


m n 


e 


h 


1.000000000 


0.9957673366 


0.9954294463 


0.9954234131 


0.9989373297 


m e 


0.9957673366 


1.000000000 


0.9996433868 


0.9996224521 


0.9904731204 


m p 


0.9954294463 


0.9996433868 


1.000000000 


0.9999732991 


0.9901455374 


m n 


0.9954234131 


0.9996224521 


0.9999732991 


1.000000000 


0.9901469965 


e 


0.9989373297 


0.9904731204 


0.9901455374 


0.9901469965 


1.000000000 



is the positive definite matrix with eigenvalues as follows: 

{4.98223, 0.0172451, 0.000495716, 0.0000263673, 6.47023 X 10 -10 }. 

Corresponding = 10 and it is close enough to our minimal accuracy, the rounding of the 
above correlator to 8 digits will make the matrix non-positive definite. 

4 As the correlation matrix of the uncertainties in the input experimental data is not a positive definite matrix 
there (supposedly by overrounding for publication), we were forced to "un-round" several matrix elements to have 
positive definite weight matrix in the least squares method of adjustment. 

6 The Plank constant is the basic one, the other selected are derived constants. In calculating the corresponding 
correlation matrix we use the minimal possible accuracy that give us the positive definite correlation matrix. 



Now we will round average values of the constants to have accuracy below the allowed 
thresholds AY . In the Table 3 we present the values of the differences (Vi) — Vf between 
calculated average values of the selected constants with the etalon accuracy and the rounded 
step-by-step values to show that after the predicted moment the end point of the rounded vector 
will be moved out of the etalon scatter ellipsoid for many standards R(V r , (V)). 



Table 3. Evolution of the "distance" of the end point of rounded vector from the etalon scatter ellipsoid 
expressed in number of standard deviations squared with rounding off the vector components 
RY in steps. 



Step 


h [J s] 


iii'e L K g| 


771 [Vcn 


"*"n L K &J 


e [CI 


R 2 (V V (V)) 


9 


4.39E-42 


2.51E-39 


1.71E-35 


4.39E-35 


-1.99E-28 


3.9E+06 


8 


3.90E-43 


-4.91E-40 


-2.25E-36 


3.91E-36 


1.4OE-30 


4.1E+04 


7 


-9.52E-45 


8.92E-42 


8.49E-38 


-8.70E-38 


1.4OE-30 


61. 


6 


4.79 E-46 


-1.08E-42 


-1.51E-38 


1.30E-38 


4.02E-31 


0.36 


5 


4.79E-46 


-7.59E-44 


4.90E-39 


3.03E-39 


2.29E-33 


0.038 


4 


-2.12E-47 


2.41E-44 


-1.01E-40 


3.16E-41 


2.29E-33 


0.00026 


3 


-1.23E-48 


4.14E-45 


-1.38E-42 


3.16E-41 


2.89E-34 


2.5E-06 


2 


-2.32E-49 


1.37E-46 


-1.38E-42 


1.62E-42 


-1.09E-35 


4.5E-09 


1 


-3.16E-50 


3.74E-47 


-3.84E-43 


-3.64E-43 


-9.03E-37 


2.7E-09 





-1.58E-51 


-2.57E-48 


1.65E-44 


3.61E-44 


9.69E-38 


6.1E-14 


AY 


45 


42 


39 


39 


31 





We see that our indices proposed as the sufficient number of digits for the safety rounding 
are indeed close to the reality. They can and should be used to the quality control of the random 
vectors obtained by statistical estimation procedures. 

Another lesson from the comparisons presented above is that the problem of the correct 
rounding off the FPC triad (Vi, Ui, Cij) is the very important problem in the task of tracing 
the space-time variability of the FPC as the improper rounding will mimic the evolution of 
constants. 

The third lesson is that the CODATA recommended values of the FPC are highly ques- 
tionable as we have convinced that the correlation matrices were corrupted by the unjustified 
rounding. 

4.2. Reliability 

As it was mentioned in the descriptions of the high level quality requirements, it is natural 
to suppose that the next iteration of the adjustment will give constants more accurate and more 
selfconsistent than the previous adjustments. 

Let us look for the time evolution of the estimates of one of the most important physical 
constant — the Planks constant h from the time of discovery up to the 2002 estimate. The 
historical perspective of the Plank constant estimates one can find in [11]. 
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Figure 1. Plank Constant: 1969-2002. Error band show that the adjustment procedures produce 
estimates that still are far from been stable, though the amplitude of variation is 
reduced in the last two releases. 



This "small-scale time variability" of the Plank constant estimates we attribute to the possible 
presence of the hidden (not estimated) systematic error introduced or missed by the adjustment 
procedures. It should be noted that systematist have to use contradictory input data which 
impossible to refine at the time of adjustment sessions 6 . 

The "evidence" of the possible stabilization (see Fig. 2 is very preliminary and should be 
tested for the other constants simultaneously by tracing the variation of the hodograph of the 
"vector of basic constants" as it is outlined in the reliability requirement. Unfortunately it is 
not possible now because of the corrupted data on correlations in the releases. The conclusion 
based on the reliability indicator is that the CODATA recommended values cannot be used in 
searches of the possible large scale space-time variations of the FPC. 



6 See discussion of this issue in the subsection: "A. Comparison of 1998 and 1986 CODATA recommended 
values" of the summary of the 1998 review ([4], pages 459-461). 
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Figure 2. Plank Constant: 1986-2002. Evidence for the possible stabilization. 



4.3. Availability 

The web access to the data on FPC offered by NIST & COD ATA in the last release (V.4.0) is 
greatly improved. Now we have easy access to all data on average values and their uncertainties 
just copy the file in the ASCII format. But unfortunately in the released list the values of 7 
basic constants out of 29 participating in the adjustment process did not quoted. The values of 
the other 28 important parameters (possible corrections to the theoretical expressions) for the 
whole adjustment procedure are omitted. They even did not discussed in the publications on 
the 1998 release. 

As it was discussed in the previous sections, the ignorance of the correlations is inadmissable 
in the high precision physics applications. But the access to the recommended correlation coeffi- 
cients remains to be the "misanthropic" one. It is hard to get data for an operative calculations 
with several constants simultaneously. There is no easy and safety way to get the complete data 
on the subsample of the triad (Vi, Ui, Cij) in a truly computer readable form. 

To extract data on say 10 constants with the correlation matrix one have to produce about 
300 flip-flops between web-pages "by hands". 

4.4. Traceability 

Traceability means that any release of the recommended FPC set should be accompanied 
with full toolkit of the input data and methods to give interested user possibility to perform all 
steps of the adjustment process and to compare the results with the recommended values. 

Unfortunately materials attached to the recommended FPC are not complete as it was 
stressed in the discussions of the availability indicator. Additional example is the incorrect 
presentation of the correlations of uncertainties in the input experimental data of the 1998 re- 
lease. 

The data on input correlations are presented only in the review on the paper [4] and the 
correlation matrix is non positive definite there [7]. 

It should be noted also that in the published documents related to the releases of FPC there 
are no discussions of the procedures used for rounding off the correlated quantities. 

5. Summary 

Summarizing the above discussions and evidences we are forced to stress that all high level 
quality requirements to the scientific information numerical data resource: correctness, selfcon- 
sistency, availability, reliability, and traceability are badly violated in at least the last three 
releases of the COD ATA recommended values of the fundamental physical constants. 

They could not be used as reference data to monitor the large scale space-time variability of 
the fundamental physical constants and moreover their usage in physics applications where the 
high precision calculations are needed is highly questionable. 

The positive outcome from our critical treatment of the quality aspects of the central numer- 
ical scientific information resource are: 
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• the preliminary proposal for the safety rounding strategy in presentation the results of high 
precision computations of the physical observables; 

• the proposal for the set of quality indicators to certify scientific information resources for 
the safety usage in physics applications; 

• the proposal of the data structure and procedures for the complete and user friendly 
Web-FPC. 
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